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Abstract. The n th row of Pascal's trinomial triangle gives coefficients 
of (1 + x + x 2 ) n . Let g(n) denote the number of such coefficients that are odd. 
We review Moshe's algorithm for evaluating asymptotics of g(n) - this involves 
computing the Lyapunov exponent for certain 2x2 random matrix products 
- and then analyze further examples with more terms and higher powers of x. 

Before discussing trinomials, let us recall well-known results for binomials. Define 
f(n) to be the number of odd coefficients in (l + x) n . Let N denote a uniform random 
integer between and n — 1, then f(N) has "typical growth" pa n 1 ^ 2 in the sense that 

E(ln(/(JV)))~iln(n) 

as n — > oo; equivalently, 

1 n ~ 1 1 
lim __£]n(/(*)) = - = 0.5. 

Also f(N) has "average growth" ps n ln ( 3 / 2 )/ ln ( 2 ) in the sense that 



ln(E(/(iV)))~^|^ln(n) 



as n — > oo; equivalently, 



n-l 

lim -^-ln I -V f(A;) I = *" v 7 ~' = 0.5849625007211561814537389. 



n^oo 



The latter value is larger since most of the Is in Pascal's binomial triangle, modulo 2, 
are concentrated in relatively few rows. Exact results are available, due to Trollope [1] 
& Delange [2] for E(ln(/(iV))) and Stein [3] & Larcher [4] for ln(E(/(A r ))); an overview 
of the subject is found in [5]. Our interest here is solely in first-order approximations. 
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Moshe [6] introduced an algorithm for evaluating such asymptotics. Define P(i,j) 
to be the (i, j) th entry of the triangle mod 2: 



/P(0,0) \ 

P(1,0) P(l,l) 

P(2,0) P(2,l) P(2,2) 

P(3,0) P(3,l) P(3,2) P(3,3) 

\P(4,0) P(4,l) P(4,2) P(4,3) P(4,4)/ 



/ 1 \ 

1 1 

1 1 
1111 
\ 1 1 / 



All values of P(i,j) in the upper right portion are Os. Define o = 0, which serves as 
a placeholder. Let £ — 1, 

AoO') = P( ,j), 0<j<£ 

and observe that the vector A = (1). Constructing additional £- vectors Ai, A 2 , . . ., 
A m _i, if required, is one aspect of the algorithm. It is mandatory that Ai(0) = 1 
always. For each < s < 1, < t < 1, let 

6^(j) = P(2o + S ,2j + t), 0<j<£ 

which is obtained by replacing (o, j) in the expression for A by (2o + s,2j + t). 
Observe that the four vectors 



(.0,0 



1,0,1 



(0), 



■i.i 



encompass only (1) = Aq and (0). No "refinement" of 6q i s hence necessary (this 
will be clarified later) and we let 

B s > t (j) = P(2o + s,2j+t), 0<j<£. 



Constructing additional £- vectors Pj'*, P^'*, 



if required, is another aspect 



of the algorithm. We have 



B, 



0,0 



Ac 



B, 



1.0 



B, 



o.i 



0. 



P, 



i,i 



Ac 



and this completes the iterative portion of the algorithm. Thus m = 1. Define mxm 
matrices via 

D s (hj) = # {i : Bf = A^ 
and observe that D = (1), D\ = (2). Define m- vectors via 



e*(j) 



1 if j = i, 
iij^i, 



w[i) 



A(j) = l} 
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and observe that eo = (1), w = (1). Let the binary expansion of a positive integer n 
be Yl r =o n r2 r with < n r < 2 and rik-i = 1. Moshe [6 J proved the following formula: 

f(n) = # {j : P(n,j) = 1} = w T D nh _ x D nk _ 2 • • • D ni D no e 

which provides a useful check that D , D 1: w are correct. 

One consequence is a well-known recursive formula for f(n). Writing binary ex- 
pansions as n = Uk-iUk-2 ■ ■ ■ nino, we see 



hence 



In = n fc _in fc _ 2 . . . nirioO, In + 1 = n k ^in k ^ 2 ■ ■ ■ nin l 

f(2n) = w T D nk _ 1 D nk _ 2 ■ ■ ■ D ni D no D e 
= i» T ViVa ' ' ' D ni D no e 
= f(n), 

/(2n + l) = w T D nk _ 1 D nk _ 2 ---D ni D no D 1 e 
= 2w T D nk _ 1 D nk _ 2 ■ ■ ■ D ni D no e 
= 2/(n). 

All _D matrices exhibited in this paper satisfy DqCq = eo, therefore even arguments 
are easy. Odd arguments are harder since .Dieo is not as predictable. 

Another consequence involves the growth rates E(ln(/(iV))) and ln(E(/(iV))). Let 
us work with the latter first. Order the complex eigenvalues fii, ■ ■ ■, of the 
matrix Dq + D\ so that \i\ has maximum modulus; therefore [6] 

ln( E (/(iV))) ; In 1^1 ln(3/2) 
ln(n) ~* ln(2) ln(2) 

as n — * oo. Working with E(ln(/(iV))) is more complicated. In this scalar case, it is 
clear that for 2 k - 1 < N < n <2 k , 

' ' " ' ' r ln(//) 



E(Hf(N))) = E [Y,HD Nr )j = Qln(l) + iln(2)) k = ^1 



ln(2) 



and so 

ln(n) ^ 2 

as n — > oo. But commutativity fails for random matrix products, in general, and 
a Lyapunov exponent-based approach will be presented in section [0.8] to deal with 
this issue. 



Odd Entries in Pascal's Trinomial Triangle 



4 



0.1. Trinomials I. Define g(n) to be the number of odd coefficients in (1 + x + 
. Properties of g(n) and Pascal's trinomial triangle are given in [7J [8]. Define 



2\n 



X 

P(i,j) to be the (i,j) th entry of the triangle mod 2: 

/ P(0,0) P(0,1) 

P(1,0) P(l,l) P(l,2) P(l,3) 

P(2,0) P(2,l) P(2,2) P(2,3) P(2,4) P(2,5) 

y P(3,0) P(3,l) P(3,2) P(3,3) P(3,4) P(3,5) P(3,6) P(3,7) ) 



\ (X \ 

1110 
10 10 10 

yiioioiioy 



All values of P{i,j) in the upper right portion are Os. Let I = 2, 

A 0') = P(o,i), 0<j<£ 

and observe that the vector A = (1, 0). For each < s < 1, < t < 1, let 

6^(i) = P(2o + S ,2i + t), 0<j<£ 

which is obtained by replacing (o, j) in the expression for A by (2o + s, 2j + t). 
Observe that the four vectors 



;i,o), 



1,0 



(1,1), 



(,0,1 



(0,0), 



:i,o) 



1,0 



contain a nonzero vector (1, 1) ^ An- Let A\ = b, 
is necessary (this will be clarified soon) and we let 

B s ' t (j)=P(2o + s,2j+t) J 0<j< 

It follows that 



;i,l). No "refinement" of b 



s.t 



B 



0,0 



Ac 



B, 



1.0 



B 



0,1 



0. 



B 



1.1 



Ar 



but we are not yet done (because of Ai). Let 

bl'Xj) = P(2(2o + s) + 1, 2(2j + 1) + 0), < j < £ 

A x by (2o + s,2j+t). 



which is obtained by replacing (o, j) in the expression for 6q° 
Observe that the four vectors 



6f° = (1, 0), 6i' u = (l,0), &r = (l,0), ^ = (0,1) 



1,0 



0.1 



1,1 



contain a nonzero vector (0,1) 7^ Ao,^. Setting A 2 = b^ 1 , however, violates the 
mandate that Ai(0) = 1 always. We thus refine the definition of b^: 

R s,t ( -, _ J P(2(2o + s) + l, 2(2j + 1) + 0) if (s, t) ^ (1, 1), 

1 [3) 1 P(2(2o + a) + 1, 2(2(7 + c) + t) + 0) if (s, t) = (1, 1) 
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where c = 1, which corresponds to shifting one step to the right. It follows that 

B®' = A , E>1' = A , B Q { 1 = A , B^ 1 = A 

and this completes the iterative portion of the algorithm. Thus m = 2. The defini- 
tions ofmxm matrices 

and m- vectors eo = (1, 0), w = (1, 2) are as before. Likewise, 

g( n ) = # {j P(ji,j) = 1} = w T D nk _ 1 D Tlk _ 2 ■ ■■D nx D nQ eQ 

as before, where n k -\n k -2 ■ ■ ■ UiUq is the binary expansion of n. 
It is easy to see that g(2n) = gin). From 

2n + 1 = n fe _in fc _ 2 • • • nin l, D x e = e + e x , 

An + 1 = n k ^ Y n k _ 2 ■ ■ ■ nin 01, D Die = 3e , 
An + 3 = n fe _in fc _ 2 • • • rcin ll, DiD^q = 3e + ei 
we reproduce Sillke's result [9] that 

g{An + 1) = 3 5 (n), <?(4n + 3) = g(2n + 1) + 2^(n). 
Also, the maximal eigenvalue fx± of Do + Di is 1 + v^5; therefore 

ln(E(o(iV))) ln(l + y/E) ln(<z>) 
v /,\ -> , , - 1 = -r4^-r = 0.6942419136306173017387902... 
ln(n) ln(2) ln(2) 

as n — > oo, where ip is the Golden mean. This constant is not new: see |10j . 

0.2. Quadrinomials. Define g${ri) to be the number of odd coefficients in (1 -+ 
x + x 2 + x 3 ) n . This extends our earlier definitions f = gi and g = g%. Define P(i,j) 
to be the (i, j) th entry of Pascal's quadrinomial triangle mod 2: 

P(0,0) P(0,1) P(0,2) \ _ / 1 

P(1,0) P(l,l) P(l,2) P(l,3) P(l,4) P(1,5)J~V1 1110 

All values of P(i,j) in the upper right portion are 0s. Let £ = 3, 

A (j) = P(o,j), 0<j<£ 
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and observe that A = (1, 0, 0). For each < s < 1, < t < 1, let 

b s ' t (j) = P(2o + s,2j + t), 0<j<£ 
and observe that the four vectors 

C = (1, 0, 0), 6j'° = (1, 1, 0), ftg- 1 = (0, 0, 0), bl' 1 = (1, 1, 0) 

contain a nonzero vector (1,1,0) ^ A . Let A 1 = b^ = (1,1,0). No refinement of 
6q* is necessary and we let 

B s ' t (j)=P(2o + s,2j+t), 0<j<l 

It follows that 

B° >° = A , Bl fi = A x , ^ = 0, B 1 > 1 = A 1 
but we are not yet done. Let 

bfU) = P(2(2o + s) + l, 2(2j +t) + 0), 0<j<£ 

which is obtained by replacing (o,j) in the expression for 6q° = A\ by (2o + s, 2j +t). 
Observe that 

C = (1, 0, 0), b{'° = (1, 0, 1), b°{ 1 = (1, 0, 0), bY = (0, 0, 0) 

contain a nonzero vector (1,0, 1) ^ A , A x . Let A 2 = b{'° = (1,0, 1). No refinement 
of b s { 1 is necessary and we let 

B'Sij) = P(2(2o + s) + l, 2(2j +t)+0), < j < t. 

It follows that 

B°'° = Ao, Bl>° = A 2 , B^^Ao, B 1 / = 

but we are not yet done. Let 

b s /(j) = P(2(2(2o + s) + 1) + 1, 2(2(2j + t) + 0) + 0), < j < I 

which is obtained by replacing (o,j) in the expression for = A 2 by (2o + s, 2j +t). 
Observe that 

bf = (1,1, 0), b 1 / = (1, 0, 1), bf = (0, 0, 0), b 1 / = (1, 0, 1) 
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encompass only Ai, A 2 and 0. No refinement of is hence necessary and we let 

B s /(j) = P(2(2(2o + s) + 1) + 1, 2(2(2j + t) + 0) + 0), < j < £. 
It follows that 



B°'° = A U Bl'° = A 2 , B°/ = 0, B? = A 



2 



and this completes the iterative portion of the algorithm. Thus m — 3. The defini- 
tions ofmxm matrices 





r 1 


2 













D = | 








1 


1. A=l 







i 




V o 





I ) 




U 


1 





and m-vectors e = (1, 0, 0), w — (1, 2, 2) are as before. Likewise, 

93(71) = # {j : P(n, j) = 1} = w T D nk _ l D nk _ 2 ■ ■ ■ D m D no e 

as before, where n k -\n k -2 ■ ■ ■ ^1^0 is the binary expansion of n. 
It is easy to see that #3(271) = g 3 (n). From 

2n+ 1 = 7ijfc_i7ijfc_ 2 . . .riiriol, £>ie = 2ei, 

4n + 1 = n k -in k -2 ■ ■ • nin 01, D Die Q = 4e , 

4n + 3 = n fc _in fe _ 2 • • • ni»oll, A-Die = 2e 2 , 

8n + 1 = n k _in k _ 2 ■ ■ • nin 001, DqDqD^q = 4e , 

8n + 3 = n fe _in fc _2 • • • nin 011, D DiDie Q = 2e±, 

8n + 5 = njt_infc_2 • • • nin 101, DiDqD^q = 8e±, 

8n + 7 = n fe _in fc _ 2 • • • nin lll, D 1 D 1 D 1 e = 4e 2 

we deduce that 

(7 3 (8n + 1) = #,(4n + 1), # 3 (8n + 3) = g 3 (2n + 1), 

(? 3 (8n + 5) = 4g 3 {2n + 1), # 3 (8n + 7) = 2g 3 {4n + 3). 

Also, the maximal eigenvalue /xi of _D + -Di is 3; it is interesting that the same 
constant 

m(EMAQ)) ; ln(3/2) 
ln(n) ~^ ln(2) 

appears here as for gi(N). 
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0.3. Trinomials II. Define h 3 (n) to be the number of odd coefficients in (1 + x + 
x 3 ) n . This extends our earlier definition g = Define P{i,j) to be the (i, j) th entry 
of the associated triangle mod 2: 

/ P(0,0) P(0,1) P(0,2) \ / 1 

VP(1,0) P(l,l) P(l,2) P(l,3) P(l,4) P(1,5)J \1 1 1 

All values of P{i,j) in the upper right portion are Os. Let £ = 3, 

A (i) = P(o,i), 0< j<£ 

and observe that A = (1, 0, 0). For each < s < 1, < t < 1, let 

bt; t (j) = P{2o + s,2j + t), 0<j<£ 

and observe that the four vectors 

C = (1, 0, 0), 6j'° = (i, 0, 0), b*' 1 = (0, 0, 0), &J' 1 = (i, i, 0) 

contain a nonzero vector (1,1,0) ^ A . Let A x = b^ 1 = (1,1,0). No refinement of 
6q * is necessary and we let 

B s ' t (j)=P(2o + s,2j+t), 0<j<t 

It follows that 

B° '° = A , B^ = Ao, P 0,1 = 0, Po 1 ' 1 ^^ 

but we are not yet done. Let 

b s /(j) = P(2(2o + s) + l, 2(2j + t) + 1), < j < t 

which is obtained by replacing (o,j) in the expression for 6q 1 = A 1 by (2o + s, 2j +t). 
Observe that 

C = (1, 0, 0), b{>° = (1, 1, 1), = (1, 0, 0), bY = (0, 1, 0) 

contain nonzero vectors (1, 1, 1), (0, 1,0)^ A , A x . Let A 2 = b{>° = (1,1,1). Setting 
A 3 = bl' 1 , however, violates the mandate that Aj(0) = 1 always. We thus refine the 
definition of b s { 1 : 

- I P ^ 2 ° + s ) + X > 2 ^ +*) + !) if («, ^ (1, 1), 

^ [J) \ P(2(2o + a) + 1, 2(2(j + c) + t) + 1) if (a, f) = (1, 1) 
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where c = 1. It follows that 

:>0,0 



Bi' 1 = A , B[ 



1.1 



An 



Br = A , B\ fi = A 2 , 

but we are not yet done. Let 

b s /(j) = P(2(2(2o + s) + 1) + 1, 2(2(2j + t) + 0) + 1), 0<j<£ 

which is obtained by replacing (o,j) in the expression for = A 2 by (2o + s, 2j +t). 
Observe that 



(.0,0 



(1,1,0), 



i,o 



(1,0,1), 6^ = (1,0,0), 6^ = (0,0,1) 



1,1 



s,t 



contain nonzero vectors (1, 0, 1), (0, 0, 1) ^ A , A u A 2 . Let A 3 = b]' = (1,0,1). 
Setting A 4 = 62' 1 , however, violates the mandate that A^O) = 1 always. We thus 
refine the definition of b s 2 

B S /(J) = 



P(2(2(2o + s) + 1) + 1, 2(2(2 j + t) + 0) + 1) if (s, t) ^ (1, 1), 

P(2(2(2o + s) + 1) + 1, 2(2(2(j +c )+t) + 0) + l) if (s, t) = (1, 1) 

where c = 2. It follows that 

B^ = A U B? = M, B°/ = A , B l / = A, 

but we are not yet done. Let 

bf{j) = P(2(2(2(2o + s) + 1) + 1) + 1, 2(2(2(2j + *) + 0) + 0) + 1), < j < £ 

which is obtained by replacing (o,j) in the expression for b 2 = A 3 by (2o + s, 2j + £). 
Observe that 

bl>° = (1,1, 0), 6j-° = (1, 1, 0), b^ 1 = (0, 0, 0), b 1 / = (1, 0, 1) 
encompass only Ai, A 3 and 0. No refinement of 63* is hence necessary and we let 
B'/ti) = ^(2(2(2(2o + s) + 1) + 1) + 1, 2(2(2(2j + t) + 0) + 0) + 1), < j < t. 
It follows that 

ti} fi = A u Bl fi = A u B°/ = 0, B\' l = A, 

and this completes the iterative portion of the algorithm. Thus m — 4. The defini- 
tions ofmxm matrices 



D = 



^12 10^ 

1 1 



\ / 



D 1 



( 1 1 1 \ 

10 1 

10 

\ 1 1 / 
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and m- vectors e = (1, 0, 0, 0), w — (1, 2, 3, 2) are as before. Likewise, 

hz{n) = # {j : Pin,]) = 1} = w T D nkl D nk2 ■ ■ ■ D ni D no e 

as before, where rik-\nk~2 ■ ■ ■ ?^i^o is the binary expansion of n. 

It is easy to see that h 3 (2n) = h 3 (n). Omitting details, we deduce that 

h 3 (An + 1) = 3h 3 (n), h 3 (8n + 3) = h 3 (2n + 1) + Ah 3 (n), 

h 3 (16n + 7) = h 3 (8n + 3) + fc 3 (2ra + 1) + 3h 3 (n), 
h 3 (16n + 15) = 2h 3 (8n + 7) + /i 3 (2n + 1) - 2h 3 (n). 
Also, the eigenvalues of D + L>i have minimal polynomial £ 4 — 3£ 3 — 2£ 2 + 2£ + 4 and 
HE(h(N))) 



ln(n) 



0.7274509132400228143266172. 



as n 



oo. 



0.4. Quintinomials. Define ^(n) to be the number of odd coefficients in (1 + 
x + x 2 + x 3 + x A ) n . Define P(i,j) to be the (i, j) th entry of Pascal's quintinomial 
triangle mod 2: 

P(0,0) P(0,1) P(0,2) P(0,3) 

P(1,0) P(l,l) P(l,2) P(l,3) P(l,4) P(l,5) P(l,6) P(1J) 

Let £ = 4 and Aq — (1, 0, 0, 0). We simply summarize: 



10 

111110 



i,o 



(1,1,1,0), A 2 = 6J> 1 = (1,1,0,0); 



5°'° - 




B i,o 


= A U 




o, 


n o 


= / 










A 3 


= b{'° = 


(1,1,1,1); 










A 


2, 




= A 3 , 


fi?' 1 = 


A), 


Bl' 1 


= A, 




A , 


ol,0 


= A , 




A , 




= A 


5°'° - 
-°3 — 


A 


2, 




= A 


2, 


^ = 


A 2 , 




= A 2 






( 1 


1 2 


\ 






(° 


1 2 


\ 




A) = 










1 



2 




A = 


1 
1 







2 








\o 





0/ 








1 







and iw = (1,3,2,4). 

The eigenvalues of Dq + D\ have minimal polynomial £ 4 — £ 3 — 6£ 2 — 4£ — 16 and 
ln(E(^ 4 (iV))) 



ln(n) 



0.7896418505307685639015472... 



as n 



OO. 
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0.5. Trinomials III. Define h^n) to be the number of odd coefficients in (1 + 
x + x 4 ) n . Define P(i,j) to be the (i,j) th entry of the associated triangle mod 2: 



P(0,0) 
P(1,0) 



P(0,1) 

f(l,l) 



P(0,2) 
P(l,2) 



P(0,3) 
P(l,3) 



P(l,4) P(l,5) P(l,6) P(1J) 

Let £ = 4 and A = (1, 0, 0, 0). We simply summarize: 



o o 
o o 



10 



D n = 







A, = C = 


(1,0,1,0); 






B°'° - 


~- A , 


B 1 / = A, 


sr = 0, 


B i r 


= A ; 


A 2 ~- 


= b°{° 


= (1,1,0,0), 


A 3 = b{'° = 


(1,1,1,1); 


B°'° = 


-A 2 , 


= A 3 , 


B[+ = 0, 


B r 


= A 2 ; 


B°'° - 
r> 2 — 


A , 


A, = b 1 / = 
B 1 / = A 4 , 


(1,1,1,0); 
B / = A , 


B r 


= A 


sr = 


A 2 , 


A 5 = 6j-° = 
B 1 / = A*, 


(1,0,0,1); 
B°/ = A 2 , 


Br 


= A 2 


B° 4 ° = 


A 2 , 


A 1.1,0 

A 6 = V = 
-B4' — Aq, 


(1,0,1,1); 
-B4' 1 — ^0, 


Br 


= A 2 


B o,o = 


A , 


A 7 = b 1 / = 
B l 5 ° = A , 


(1,1,0,1); 
B^ 1 = A , 


sr 


= A 7 


B°'° - 
D % — 


Ai, 


5a 1 ' = A 7 , 


Bq A = Aq, 


sr 


= A 5 


B° 7 >° = 


A , 


B 1 / = A 2 , 


Br = a 2 , 


B i r 


= A A 



( 


1 





2 





1 


2 


1 


1 


\ 



































1 





2 


1 





1 


1 














































































































































/ 



( 


1 





1 








1 












1 






























1 





1 


1 








1 









1 






























1 














1 















1 








1 





















1 













V 

















1 


1 





) 



and w = (1,2,2,4,3,2,3,3). 

The eigenvalues of Dq + D\ have minimal polynomial £ 5 — 3£ 4 — 2£ 2 — 8£ + 8 and 

ln(E(/i 4 (N))) 



ln(n) 



0.7362115557393079316549209. 



as n 



00. 
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0.6. Sextinomials. Define g^(n) to be the number of odd coefficients in (1 + x + 
. . . + x 4 + x 5 ) n . Define P(i,j) to be the (i,j) th entry of Pascal's sextinomial triangle 
mod 2: 



P(0,0) P(0,1) 
P(1,0) P(l,l) 



P(0,4) 

P(l,4) P(l,5) P(l,6) 



P(l,9) 



1 

1111110000 



Let £ = 5 and A = (1, 0, 0, 0, 0). We simply summarize: 



A 1 



1,0 



10, 



A 1,0,0 
A 2 = V = 



(1,1,1,0,0); 



^o ' 1 = 0, 



B, 



1,1 



^i 



0,0 



5 2 °'° = A 



0, 



(1,1,0,0,0); 
B\'° = Ai, B®' 1 = A , 
A 3 = b 1 / = (1,0, 0,1,0); 

i^ = A 3 , flj 1 = A3, 



si- 1 = a i; 



So 1 ' 1 = 



5. 



0,0 



A, 



0, 



A 4 
B 



1,0 



1,0 



A, 



A 5 = 6°'° = (1,0, 1,0,0); 



(1,1,0,1,1); 



B 



o.i 



An 



5 



BT = A, 
B° 5 ° = A 2 , 



Bi° = As, 



B / = A 2 , 



B 



1,0 



^0 



/ 1 1 2 2 \ 



10 11 







Bl l = 0, 



D 1 



B 
B 



i.i 



1,1 



0; 

As; 



4 

1.1 



A 3 ; 

5 

/000000\ 
2 2 

112 2 
1 
\000000/ 



\ 1 / 

and w = (1,3,2,2,4,2). 

The eigenvalues of D +Di have minimal polynomial £ 6 — 4£ 5 +£ 4 — £ 3 + 8£ 2 + ll£+8 
and 

HeMn))) 



ln(n) 



0.8194694621655401465959376... 



as n 



oo. 
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0.7. Septinomials. Define ge(n) to be the number of odd coefficients in (1 + x + 

hi 5 + x 6 ) n . Define P(i,j) to be the (i,j) th entry of Pascal's septinomial triangle 

mod 2: 



P(0,0) P(0,1) 
P(LO) P(l,l) 



P(0,5) 

P(l,5) P(l,6) P(l,7) 



P(l,ll) 



1 

111111100000 



Let £ = 6 and A = (1, 0, 0, 0, 0, 0). We simply summarize: 



A 1 



(1,1,1,1,0,0), A 2 = 6 ' 1 = (l,l,l,0,0,0); 



R°'° - A 



Br = ^ 



.1,0 _ A 
— -^1) 



(0,0 



- ^3, 



- D 2 



5; 
5 



0,0 



As, 
Ao, 



A 3 = b» 

Br 

a, = b i r 

A 4 , 
Ao, 



B 



o,i 



0, B t 



1,1 



(1,1,0,0,0,0); 



B°r = A 3 , 



= A 2 ; 
B i r = A 3 ; 



br = (1,1, 1,1, 1,0); 



□ 1,0 



B. 



1,0 



B 
B 



o,i 

2 

0,1 
3 



Ao, 
Ao, 



B. 
B. 



Li 

2 

1,1 



A 3 ; 
Ao; 



A, 



1,0 



(1,1,1,1,1,1); 



0,0 

4 — 


A 2 , 


Br 


= A 


5, 


B°r = 


As, 




B i r 




Ao 




0,0 

5 ~~ 


A 2 , 


sr 


= A 


2, 


B°r = 


A 2 , 




Bir 




A 2 


> 


/ 1 


1 


2 


\ 






( 








2 


1 


\ 


















1 























1 


2 




D x = 


1 














2 





2 1 


1 










2 


1 

































1 











\o 








oy 






\o 











1 


0/ 



Dn 



and w = (1,4,3,2,5,6). 

The eigenvalues of D + Di have minimal polynomial £ 6 — £ 5 — 2£ 4 — 28£ 3 
and 

ln(E(<? 6 (A0)) 



ln(n) 



0.8317963967344406899938931... 



16^ + 64 



as n 



oo. 
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0.8. Lyapunov Exponents. All Dq matrices shown in this paper satisfy rank(/}p) 
1 for some positive integer q; further, there is a change of coordinates under which Dq 
is transformed into the matrix whose (0, 0) th entry is 1 and all of whose other entries 
are 0. That is, there is an invertible m x m matrix Q with Q~ 1 D$Q = e§e^. Define 

D' Q = Q^DoQ, D[ = Q- 1 D 1 Q 

and let z , zi, . . ., z k _ 2 , z k _i denote a sequence of independent random coin tosses 
(heads=l and tails=0 with equal probability). The Lyapunov exponent corresponding 
to random products of D' and D[: 



A = lim — In 

fc— >oo k 



D' D' ■ ■ ■ D' D' 



exists almost surely. It turns out that A/ ln(2) is precisely what we seek to characterize 
typical growth rates of gi{N) and hj(N). 

Let x(0 9 ) denote the set of all finite binary words z with no subwords q and with 
rightmost digit 1. (0 1 means 0; 2 means 00; 3 means 000.) Let £(z) denote the length 
of z. For example, 1111 is the only word of length 4 in %(0); 0101, 0111, 1011, 1101, 
1111 are the only words of length 4 in x(00)> the set x(000) additionally contains 0011 
and 1001. It is natural to sort the elements of x(0 9 ) m terms of increasing length. 
Write z = z oZl . . . z k _ 2 z k _ x and D' z = D' zo D' zi ■ ■ ■ D' Zk _D' Zk x . By D' z (0, 0) is meant 
the upper left corner entry of D' z . 

Extending earlier work by Pincus [11] and Lima Sz Rahibe [12] , Moshe [13] proved 
that 

A = — Try- r V -^r In \D' (0, 0) I . 

v ; ze x (o q ) 

This series is attractive, but computationally difficult since the number of words in 
x(0 9 ) of length k grows exponentially with increasing k. Summation of the series, 
coupled with Wynn's ^-process for accelerating convergence, serves as our primary 
method for calculating A. Our secondary method is based on the cycle expansion 
method applied to a corresponding Ruelle dynamical zeta function [TH [151 EE EZ] • 
In the case of the first trinomial (1 + x + x 2 ) n , we have q = 1, 

Moo)' ^{\~-\ 

hence [6] 

k=l k=l v 7 

= 0.4299474333424527201146970... 
= ln(1.5371767171823579495901403...) 
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and 



E(Hg 2 (N))) 



A 



ln(n) ln(2) 
asn-> oo. This is smaller than \n(ip)/ ln(2) 



0.6202830299260946960737425... 



In the case of the quadrinomial 



Q 



= 0.694..., as discussed earlier. 
1 + x + x 2 + x 3 ) n , we have q = 2, 



1 


-2 






( 1 















1 


|, 1*0=1 











1. «-( 





1 


I J 




V o 


1 








-6 



hence 



x = - V JL 

24 ^ 2 £ W 

zex(oo) 



In |^(0, 0)| = 0.34657359... 



and A/ln(2) = 0.49999999... < 0.584.... We conjecture that A/ln(2) equals 1/2 (the 
typical growth rate for binomials) and prove this to be true in section |0.12j . The 
cycle expansion converges slowly in this case, therefore Moshe's technique is helpful 
here. 

In the case of the second trinomial (1 + x + x 3 ) n , we have q — 2, 

/ 1 -3 -2 1 \ ^10 






Q 





1 





1 










D> 



J 





1 








0/ 



/3 


1 



-6 






hence A = 0.45454538229305... and A/ln(2) 
related example appears in [13] . 

In the case of the quintinomial (1 + x + x 2 + x 3 + x 4 ) n , we have q 

( 1 -3 -2 2 \ /1000\ 



1 



0.65577036889316.. 



1 
< 0.727... 



4 \ 

2 

0/ 

A 



01 -2 0000 
00 1 ' °~ 0100 
\00 1/ \0000/ 

hence A = 0.504253705692... and A/ln(2) = 0.727484320552... < 0.789.... 



/ 5 
1 
1 



-10 

-1 

-3 
1 



-8 4 
-2 -2 
-2 4 



In the case of the third trinomial (1 + x 



x 



4\n 



we have q = 2, 



/ 


/ 


1 


-2 


-2 





-1 


-2 


-1 


-1 


\ 


\ 









1 





-2 


-1 





-1 


-1 
















1 


































1 


































1 


































1 


































1 









V 


V 























1 


/ 


/ 
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D' 



( 


1 








o 


o 


o 


o 


o 






( 


3 





-2 


-8 


-4 


-2 


-4 


-4 


\ 


































1 





-1 


-4 


-2 


-1 


-2 


-2 









1 





























1 





-1 








-1 










































1 





-2 


-1 





-1 


-1 




















































































1 














1 













































1 








1 



















































1 













\ 


























J 




V 

















1 


1 





/ 



hence A = 0.45759385431410... and A/ln(2) = 0.66016838436022... < 0.736.... 



In the case of the sextinomial (1 



x 



Q = 



hence 





1 


-1 


-2 


-2 


-1 


-2 






( 


i 

















\ 







-1 








1 










































1 

















1 
































, D' = 



































1 














1 


















1 








































V 








1 





-1 





J 




V 




















/ 



x + x ) n , we have q 



3, 

( 6 -8 



2 -4 



2 -4 

y o 2 




-4 

-4 

2 



-10 

-4 

-4 
1 



-4 





-6\ 
1 

-3 








-2 



A 



1 

112 



1 



ln|D;(0,0)| = 0.5344481528... 



zex(ooo) 

and A/ln(2) = 0.7710456996... < 0.819.... 

In the case of the septinomial (1 + x + • • • + x 5 + x 6 ) n , we have q 



3, 



Q = 



( 


1 


-6 


-6 


-4 


1 





\ 



















1 












2 





-1 


















2 
























1 


-2 




\ 





1 








1 

"2 


1 


/ 



D' = 



/iooooo\ 



1 

1 



yooooooy 



f 


7 


-36 


-28 


-24 


4 


-4 


\ 










1 





1 
2 


-2 






3 
2 


-8 


-8 


-6 


1 

2 
1 

"2 


1 












1 





1 






2 


-12 


-10 


-8 


1 







\ 


1 


-6 


-6 


-4 


1 





/ 



hence A = 0.53765282... and A/ln(2) = 0.77566905... < 0.831.... 

Let us conclude by mentioning relevant sequences in Sloane's online encyclope- 
dia. A001316 is f{n); A000120 is ln(/(n))/ln(2); A006046 is £ fc<n /(£;); A000788 
is £*<„M/(*0)/M2). A071053 is g(n); A134659 is E fc <„^)- A134660 is g 3 (n); 
A036555 is la(g 3 (n))/ ln(2). A134661 is h 3 (n). A134662 is h 4 (n). A007318, A027907, 
A008287, A035343, A063260, A063265 are Pascal's triangles associated with (l + x + 
■■■ + x 1 " 1 + x r ) n for r = 1, . . . , 6; A038717 and A134663 are likewise for (1 + x + x 3 ) n 



and (1 



. It is well-known that ln(/(n))/ln(2) is the number of Is in the 



binary expansion of n, but scarcely noticed that ln(g 3 (n))/ ln(2) is the number of Is 
in the binary expansion of 3n. 
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0.9. Pascal's Rhombus. The sequence of polynomials giving Pascal's trinomial 
triangle arises from the first-order recurrence 

p n (x) = (1 + x + x 2 )p n _i(x), p (x) = 1. 

Pascal's rhombus [18], by contrast, arises from the second-order recurrence 

p n (x) = (1 + x + x 2 )p n ^i(x) + x 2 p n _ 2 (x), px{x) = 1 +x + x 2 , p (x) = l. 

In this addendum, we perform the same analysis as in the preceding. Define u(n) to 
be the number of odd coefficients in p n (x). Define P(i,j) to be the (i, j) th entry of 
the associated "rhombus" mod 2: 



P(0,0) P(0,1) P(0,2) P(0,3) 

P(1,0) P(l,l) P(l,2) P(l,3) P(l,4) P(l,5) 

P(2,0) P(2,l) P(2,2) P(2,3) P(2,4) P(2,5) P(2,6) P(2,7) 



10 
1110 
10001000 



The objects A± and B t are no longer ^-vectors but 2 x i matrices. It is mandatory 
that Ai(0, 0) = 1 or 0) = 1 (or both) for every i. Let t = 4 and 



We simply summarize: 



At 



t.0,0 



An 



10 
10 10 



10 
1110 



A, 



B 



0,0 



Ax, 



B, 



1,0 



A, 



As = K 



0,1 



Bl 



o.o 




110 

>i,o _ 



B 
B 



o.o 



0.0 



B 



0.0 



A), = A 2 , 

Ax, Bl'° = A , 

A,, Bl>° = A,, 

A 3 , Bl'° = A, 

( 1 \ 

1 2 


10 1 

\ 2 1 ) 



• i,o 



110 
10 1 



o.i 



B, 



A 4 = Bl' 1 



B^ 1 = A 3 , Bl 



B, 



1,1 



A ; 




10 



i,i 



B^ 1 = A,, 



B 



2 

0,1 



A 



■i- 



B. 
B. 



i.i 

2 

1,1 



B 



0,1 



Aa 



B 



i.i 



A 4 ; 
A ; 
A\, 
= 0; 



/ 1 2 \ 

2 1 

110 



\ 1 / 
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and w = (1, 1, 2, 0, 0). The maximal eigenvalue of D + Di is (3 + VT7)/2; therefore 

MEKAO)) ^ ln((3 + 717)74) = . 8 3 2506 3 8 3 5804514437981667 ... 
ln(n) ln(2) 

as n — > oo. To compute lim^oo E(ln(w(iV)))/ ln(n) via Moshe's technique requires 
a binary word 2 for which the product D z satisfies rank(D z ) = 1. No such word 
z exists, therefore our primary method is inapplicable here. Our secondary method 
gives A = 0.57331379313..., hence A/ln(2) = 0.82711696622.... 
We mention the Fibonacci polynomials 

p n (x) = xp n „i(x) +p n - 2 (x), Pi(x) = x, p (x) = 1 

and that the number v(n) of odd coefficients in p n (x) is the n th term of Stern's 
sequence [TH [201 HI] 

v{2n + 1) = v(n), v{2n) = v{n) + v{n — 1). 

It follows that 

Do= (\ ?)' Dl = (li)> « , = ( 1 >°); 

ln(E(«(iV))) ln(3/2) E(ln(«(JV))) A 

,\ 1,1 -> \\ , A\ -> r— - = 0.571613901650254... 

ln(n) ln(2) ' ln(n) ln(2) 

as — >• oo. The estimate A = 0.396212564297744... is new (as far as we know), but 
the challenge of computing A was posed long ago [22] . 

Relevant sequences in Sloane's online encyclopedia are A059319 for u(n), A059317 
for Pascal's rhombus, A002487 for v(n) and A049310 for "Fibonacci's rhombus". We 
recall Glaisher's theorem that f(n) is the number of odd binomial coefficients of the 
form ("J, < m < n; this is mirrored by Carlitz's theorem that v(n) is the number 
of odd binomial coefficients of the form ( n ~ m ) Q < 2m < n. 

0.10. Extreme Values. The function f(n) has "maximum growth" pa n 1 in the 
sense that 

limsup — - In(/(n)) = 1. 
To evaluate the limit superior, we note that 



f(n) < f(2 k -l) = 2 k for all n < 2 k - 1 
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for each k = 0, 1, 2, . . .. No other disjoint subsequence possesses this property and 

ln(2 k ) 

iim - — 7-7 r = 1. 

fc->oo ln(2 fc - 1) 
The functions g(n) and u(n) likewise satisfy 

1 



limsup 



ln(n^ 

In the case of g(n), we have 



ln(g(n)) = limsup 



\n(n 



■ ln(w(n)) 



g(n) < 



9C2 k - 1, 
0(3 • 2 k+1 - 1 



2 k+2 - (-1) 



for all n < 2 k — 1, 
for all n < 3 • 2 fc+1 - 1, 



#(11 • 2 k+l - 1) = 3(2 fe+3 + (-l) fc ) for all n < 11 • 2 fc+1 - 1 



for each fc = 0, 1, 2, . . .. No other disjoint subsequence possesses this property and 



In 



lim 



2 fc+2_(_ 1 )fc 

3 



lim 



ln(2 fe+3 + (-1) 



ln(3(2 fc+3 + (-l) fc )) = 



*-« ln(2 fc -l) fc->oo ln(3 • 2 k+1 - 1) fc^oo ln(ll • 2 k+1 - 1 
In the case of u(n), we have 



2 k 



u(n) < < 



' u(5) = 6 
u (37) = 45 

w(2 fc - 1) = 
u(5 • 2 fe+1 - 1) 
u(2 • 4 fe + 2 - 7) 



5 . 2 k+3 + (-l) k 
3 

5 • 4 fc+2 + 12k + 1 



for all n < 5, 
for all n < 37, 

for all n < 2 k — 1, 

for all n < 5 • 2 fc+1 - 1, 

for all n < 2 • 4 fe + 2 - 7 



for each = 0, 1, 2, No other disjoint subsequence possesses this property and 



In 



lim 



2 fc+2-(-l)fc 

3 



hi 



lim 



5-2 fc+3 +(-l)' 
3 



hi 



lim 



5-4 fc+2 +12fc+l 
3 



ln(2 fe - 1) fc^oo ln(5 • 2 k+1 - 1) ln(2 • 4 fc + 2 - 7) 

By contrast, the function v(n) satisfies 
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where ip = (1 + \/5)/2. Let ao = 0, ai = 1, % = a 3 -_i + 0,-2 denote the Fibonacci 
sequence. We have 



v(n) < 



2(4 fe -l)\ _ „ 2(4 fc -l) 

~ a 2k+i tor all n < 



3 J " 3 

4(4 fc -l)\ 4(4* -I] 



fl2fc+2 for all n < 



3 

for each k — 0, 1, 2, No other disjoint subsequence possesses this property and 

In(o 2 fc+i) ,. ln(a 2k+2 ) ln(<p) 

lim 7 : r- = iim 7 ; r- = -— . 

fc ^°°m ( 2{ M fc -"°° In ( 4( } ) ln ( 2 ) 

Our analyses in this addendum are aided by the -Do, D\, w matrices from earlier 
sections. "Minimum growth", defined with limit superior replaced by limit inferior, 
is not as interesting for three of the cases since 

f(2 k ) = 2, g(2 k ) = 3, v(2 k - 1) = 1 

always. The remaining case, u(n), resists all attempts at simplification. 

0.11. Variability. We shall be very brief here. Let N denote a uniform random 
integer between and n — 1. Kirschenhofer [23] proved that f(N) has "typical 
dispersion" w n ln( ^ 2 ^ 4 in the sense that 



as n — > oo; equivalently, 



lim 



ln(n) 



Var(ln(/(iV))) ~ ^2^1n(n) 



n-l n-1 \ • 

-^ln(/(A;)) 2 - 

A;=0 \ fc=0 / 



ln(2) 



In fact, this is related to a generalized Lyapunov exponent of order two corresponding 
to random products of D = (1) and D 1 = (2). We defer further study of such 
quantities to a later paper. 

Define ip — (5 + vT7)/4 for convenience. As another example, v(N) has "average 
dispersion" m 7i lxi w/ ln ( 2 ) in the sense that 

ln(Var(,;(iV)))~^|ln(n) 
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as n — > oo; equivalently, 

1 



lim 

n-+oo m{n) 



In 



1 n— 1 / n— 1 



fc=0 



fc=0 



ln(V') 



More details on the latter result can be found in [2H [251 [26l [271 [28l [29l [30] . In fact, ip 
is the leading eigenvalue of |(A) <S> A) + where ® is the direct or Kronecker 

product of matrices 



1 
1 1 



thus 



D ®D C 



( 1 \ 
110 
10 10 

Vi i i 1/ 



D l ®D 1 



I)' 








/I 


1 


1 


1 \ 





1 





1 








1 


1 










1 / 



In the same way, 

ln(Var(/(iV))) ln(5/2) 



\n{Vav{g{N))) ln(0 



ln(n) 



ln(2) 



ln(n) ln(2) ' 

where £ = 2.813... has minimal polynomial £ 3 — 2£ 2 — 3£ + 2 and 

ln(Vax(u(JV))) Hv) 
ln(n) " ln(2) 

where rj = 3.194... has minimal polynomial 4rj 7 — 8r] 6 — 25i] 5 + 22r] 4: + 24:rj 3 + lQr]' 2 +ri — 2. 

0.12. Proof of Conjecture. Since the Lyapunov exponent is defined as an almost- 
sure limit: 



A 



lim - In \\D Zo D Zl 



D Zk . 2 D Zk A 



k— >oo k 

we may assume that z = 1. Every binary word z thus looks like 

1 jo 1 jl 1 J2 ... 1 jn ~ 2 1 O^"- 1 
where each j$ > 0. We introduce a rewording of D z : 
D j0 D h D h ■ ■■!),, = (DM ) {DiD^) (D 1 D^) • • • (d 1 D 



Jn-2 




Jn-l 
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and note that the probability associated with Dj i is l/2 J ' i+1 . Clearly 



A = lim — In 

n— >oo 2n 



D j0 D h 



Da -D* , 

Jn — 2 Jn — 1 



since the mean length of Dj. is 1 + Ylj>oj/^ +1 - Also 
D„ Dj 12 I , D x = D,D 





D, 




D x Dl =[244 




for all j > 2. Hence the rewording actually consists of only three matrices Dq, D±, 
D2 with probabilities 1/2, 1/4, 1/4. Further, the initial row of each matrix is zero, 
thus we may consider only the lower-right 2x2 submatrix: 





1 2 



4 
1 



M, 



( 4 






\ o 


S H 





Let pi = 1/2, p 2 = 1/4 denote the weights corresponding to aij3f, a.<i($i and q — 1/4 
denote the weight corresponding to M. The case of two 2x2 matrices, one with 
rank 1 and the other with rank 2, was solved in [TTJ [T2]. Our case involves three 
matrices, two with rank 1 and one with rank 2, as well as the scaling factor 1/2 due 
to rewording. Generalizing, we obtain 



A 



1 2 2 oo 



=1 3=1 fc=0 



Observe that 




if i 
if z 
if z 
if i 



1, J 
1, J 



2, J 
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and therefore 

2A 1^1 1 ^ 2 1 ^ 2k 1 ~2(ife + l) 

ln(2) ~~ 4^4^ + 8^4^ + 8^ + 16^ 4 1 

v y fc=0 fc=0 fc=0 fc=0 



= 1, 



as was to be shown. An independent proof of the conjecture was found by Thomas 
Doumenc. 

Let #(n) denote the number of Is in the binary expansion of n. We know that 
#(n) = ln(/(ra))/ln(2) and 



n ! 



fc=0 

as n — > oo. A corollary of our proof is that 

Ln/3j 

21n(2) 3 



fc=0 

because #(3n) = ln(^ 3 (n))/ln(2). The formulas 

Ln/3J L«/3J 

S #(3A; + 1} ~ 2lnT2li lnH ~ S #( 3A; + 2 ) 

fc=0 ^ ' k=0 

follow similarly, that is, counting binary Is is (on average) independent of ternary 
residue. We wonder whether simpler proofs of this fact can be found. 
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